Nonequilibrium Transport in Superconductor/Ferromagnet/Superconductor Diffusive 
Junctions: Interplay between Proximity Effect and Ferromagnetism. 
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The theory of the I-V characteristics in diffusive superconductor/ weak ferromag- 
net /superconductor (SFS) junction is developed. We show that the exchange field h of the 
ferromagnet manifests itself as an additional conductance peak at e V ~ A + h in the phase-coherent 
regime, when the Thouless energy is of the order of superconducting order parameter. The excess 
current exhibits non-monotonous dependence on the exchange field and non-trivial temperature 
behavior, which is strongly influenced by the temperature dependence of the exchange field. 

PACS numbers: 74.45.+C, 74.50.+r 
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The recent progress in the experimental techniques has 
made possible the fabrication of mesoscopic structures 
on the nanometer scale. Hybrid structures containing 
superconducting and ferromagnetic elements offer an op- 
portunity to generate and control coherent spin trans- 
port with otherwise conventional electronics. The equi- 
librium transport and proximity effect in such structures 
have been theoretically and experimentally investigated 
recently in details as for the case of weak ferromagnetic 
alloys so as for half-metals like CrO^i^.In particular, 
Josephson current in SFS junctions and T c of SF bi- 
layers and multilayers have been investigated in details 
for the case of weak ferromagnetic alloys (see Ref|H and 
references therein). Equilibrium density of states was 
also studied 4 ^. On the other hand, to the best of our 
knowledge, the nonequilibrium transport in SFS meso- 
scopic junctions has not been studied yet neither from 
theoretical nor from experimental point of view except 
for theoretical investigations of magnetic quantum point 
contacts^ 7 -^. We address our theoretical paper to the 
part of this problem and study the phase-coherent trans- 
port in diffusive voltage-biased SFS plane junctions. We 
consider interlayers made of a weak ferromagnetic mate- 
rial, in which the value of the exchange field h (measured 
in the energy units) is of order of superconducting order 
parameter A. In ferromagnetic alloys like CuNi, which 
have been intensively used by now for experimental in- 
vestigation of equilibrium properties of SFS heterostruc- 
tures, the exchange field is several times larger than A. 
On the other hand, it can be concluded from our analysis 
of the problem that for studying the TV characteristics 
of SFS junctions, made on the basis of a weak ferromag- 
netic alloy, the most interesting case is h < A. As far as 
we know, the work on the creation of appropriate alloys 
is in progress now, so we believe that this limit can be 
experimentally realized in the nearest future. 

The TV characteristics of the superconductor/normal 
metal/superconductor (SNS) voltage-biased diffusive 
junction are studied theoretically in details mainly in two 
limits. The first one is the limit of short junction d -C £, 
where d is the length of the interlayer, £ = y/KD j A is the 
superconducting coherence length and D is the diffusion 
constant. In this regime the subharmonic gap structure 



(SGS) in the differential conductance dl/dV consists of 
a set of pronounced maxima at eV n — 2A/ro 1 °'- 11 . The 
excess current, which takes place at high biases eV ^> A, 
behaves in dependence on temperature as A(T) does and 
can be as positive so as negative depending on the trans- 
parency of SN interface o 10 ' 12 . The second one is the in- 
coherent limit d > ^. The proximity effect is negligi- 
ble in this regime and the transport is determined by 
the kinetic equation for the distribution functio n 13 : 14 ' 15 . 
The SGS of the conductance again shows sharp features 
at eV n = 2A/n. At the same time in the intermediate 
regime £ ~ d the proximity effect, which manifests it- 
self as a minigap A g in the equilibrium density of states 
in the normal region, takes place. It has been shown 
very recently^ that for SNS voltage-biased diffusive junc- 
tion with highly-transparent NS interfaces in this regime, 
when the interplay between proximity effect and MARs 
takes place, the well-known subgap conductance struc- 
ture eV n = 2A/n modifies and exhibits an additional 
maximum at roughly eV ~ A + A g . These predictions 
are in good agreement with the existing experiment s 17 : 18 . 
In the present work we numerically calculated the I-V 
characteristics taking into account the proximity effect 
and non-equilibrium distribution function in the ferro- 
magnetic interlayer. We show that in the phase-coherent 
regime, when the proximity effect in the ferromagnetic 
region takes place, the exchange field of a weak ferro- 
magnet h < A explicitly manifests itself in the SGS of 
diffusive voltage-biased SFS junction leading to the ad- 
ditional maximum at eV ~ (A + h), which is the most 
pronounced in the case of low-transparent SF boundaries. 
For the case h > A the differential conductance exhibits 
no pronounced characteristic features. The dc current 
at high voltage biases is also investigated. It is deficit 
(Idef = I — V/R < 0) for the case of highly- resistive in- 
terfaces we consider and has a non-monotonous behavior 
as a function of exchange field, reaching the maximum 
value at roughly h ~ A. This non- monotonous behav- 
ior results in the fact, that the temperature dependence 
of the deficit current exhibits the characteristic features, 
which are determined by the temperature dependence of 
the exchange field and can be used for its experimental 
mapping. 
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Further the model under consideration and the method 
we use are briefly described. We study an SFS junc- 
tion, where F is a diffusive weak ferromagnet of length 
d coupled to two identical superconducting (S) reser- 
voirs. The superconductors are supposed to be diffusive 
and have s-wave pairing. We assume the SF interfaces 
to be not fully transparent and suppose that the resis- 
tance of the SF boundary R g dominates the resistance 
of the ferromagnetic interlayer Ftp- We use the qua- 
siclassical theory of superconductivity for diffusive sys- 
tems in terms of time-dependent Usadel equations^. The 



fundamental quantity for diffusive transport is the mo- 
mentum average of the quasiclassical Green's function 
g(x,e,t) = {g(pf,x,E,t)) Pf . It is a 8 x 8 matrix form 
in the product space of Keldysh, particle-hole and spin 
variables. Here x - is the coordinate measured along the 
normal to the junction, t stands for a time variable and 
e is the excitation energy. 

The electric current should be calculated via Keldysh 
part of the quasiclassical Green's function. For the plane 
diffusive junction the corresponding expression reads as 
follows 
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where e is the electron charge and h — 1 throughout the 
paper. (g(x, e, t) ® (dg(x, e, t))/(dx)) K is a 4 x 4 Keldysh 
part of the corresponding combination of full Green's 
functions. The product ® of two functions of energy 
and time is defined by the noncommutative convolution 
A®B = e l[d e 9 t ~9t9f) A{e,t)B{e,t). f t and a l are Pauli 
matrices in particle-hole and spin spaces respectively. 

The quasiclassical Green's function g(x, e, t) satisfies 
the non-stationary Usadel equation. In order to solve 
the Usadel equation it is convenient to express quasi- 
classical Green's function g in terms of Riccati coher- 
ence functions -y R < A and ; y R ' A and distribution functions 
x K and x , All these functions are 2x2 matrices 
in spin space and depend on (x, s, t). The correspond- 
ing expression for g can be found in Ref[23| and there- 
fore is not explicitly written here. Riccati coherence 
and distribution functions obey Riccati-type transport 
equations 2 ^, where the appropriate self energy takes the 



form tl(x,e,t) = h&s + {l/2i:T s )cfj,gcf^. Here t^ 1 is an 
inverse magnetic scattering time. As the F layer is sup- 
posed to be an alloy, a role of magnetic scattering may 
be quite import an t 20 i 21 . We assume a presence of the 
relatively strong uniaxial magnetic anisotropy which pre- 
vents mixing of spin-up and spin-down Green's functions, 
so the magnetic scattering term is a diagonal matrix in 
spin spaced. 

The Riccati-type transport equations should be solved 
together with the boundary conditions at SF interfaces. 
As it was mentioned above we consider the case when 
the dimensionless conductance of the boundary G = 
Rp/Rg < 1, so the interface transparency T ~ G(l/d) <C 
1 (I is the mean free path). Due to the smallness of the in- 
terface transparency T we can use Kupriyanov-Lukichev 
boundary conditions at SF boundaries 25 . In terms of Ric- 
cati coherence and distribution functions they take the 
form 
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Here Riccati coherence and distribution functions de- 
noted by the lower case symbols /, r are taken at the left 
and right ends of the ferromagnet. The quantities de- 
noted by the lower case symbols (S; I, r) are correspond- 
ing Green's functions at the superconducting side of the 
left and right SF interfaces. We assume the parameter 
{Rf / Rg){<JF / &s), where up and a s stand for conductiv- 



ities of ferromagnetic and superconducting materials re- 
spectively, to be also small, what allows us to neglect the 
suppression of the superconducting order parameter in 
the S leads near the interface and take the Green's func- 
tions at the superconducting side of the boundaries to 
be equal to their bulk values with the appropriate shift 
of the quasiparticle energy due to the applied voltage 
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V = V r -Vi. 

Let us now turn to the analysis of the I-V characteris- 
tics. 
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FIG. 1: The zero temperature differential conductance 
dl/dV, normalized to its value for the normal state of the 
leads, as a function of eV/T c , taken for: (a) different ex- 
change fields: h = 0(gray solid); h — 0.45T C (dotted); h = 
0.9T C (black solid) and h = 2. 6T C (dashed); The insert to 
(a) shows the Rdl/dV as a function of eV/T c for h > A: 
h = 2.6T c (solid); h = 3T c (dashed) and h = 3. 5T C (dotted) 
are practically undistinguishable; (b) different values of the 
Thouless energy: eth = 5T C (dashed); eth = 2T C (black solid) 
and Eth = 0.5T C (dotted). T c is the superconducting criti- 
cal temperature. G = 0.1 and l/r s = 0. The inset to (b) 
demonstrates the influence of the magnetic scattering on the 
conductance: l/r s =0 (solid); 1/t s — 0.05T C (dotted); l/r s = 
0.15T C (dashed); l/r s = 0.3T C (dashed-dotted); l/r s = 0.45T C 
(gray solid); l/r s = 0.9T C (gray dashed), h = 0.9T C . 

Figure QJa) shows the dependence of the zero tempera- 
ture differential conductance dl/dV on voltage for differ- 
ent values of the exchange field in the ferromagnet. The 
following main features can be observed: 

(i) there is a series of peaks at eV n — 2A/n. In fact, 
only the peaks corresponding to n = 1, 2 are seen in the 
figure. The other harmonics are smeared out due to the 
small interface transparency and large enough intrinsic 
broadening of the quasiparticle energy (modelling the in- 
elastic scattering rate) 7 = 0.03T C , which we take in the 
numerical calculations. This value of intrinsic broadening 
is large, but consistent with an estimate of inelastic scat- 
tering rate due to electron-phonon processes 7 ~ e /u^ 
at high enough energies e ~ 



(ii) for h < A the large additional maximum appears 
at roughly eV — A + h, while the curves corresponding to 
h > A do not exhibit any additional pronounced charac- 
teristic features, in particular at eV = A + h (See Insert 
to FigfJJa)). The exchange fields of order of several A 
only result in the suppression of the peak at eV = A. 
In fact, the larger the exchange field the less different 
from each other the appropriate conductance curves. Of 
course, this statement is related to the region of weak 
exchange fields we study. 

Analyzing the dependence of the I-V characteristics on 
the Thouless energy etk = D/cP, which is represented in 
Figure QJb), we can conclude, that the most pronounced 
maximum at eV ~ A + h exists at the region Eth ~ 
A. If the Thouless energy increases the peak washes out 
and finally vanishes. At the limit of very short junction 
Erh 3> A the only structure eV — 2A/n survives. On 
the other hand, with decreasing the Thouless energy the 
amplitude of maximum at eV ~ A + h reduces to the zero 
value at the limit of very small Thouless energies and 
again the only structure at eV = 2A/n can be observed. 
With the increasing of transparency of the interfaces the 
maximum smears out and is masked by the growing peak 
at eV = A. 

The inset to Figure QJb) demonstrates the influence of 
the magnetic scattering on the maxima discussed above. 
It can be seen that the magnetic scattering does not in- 
fluence the peak at eV = 2A, while the maximum at 
eV ~ A + h is broadened, but not dramatically for small 
compared with T c values of magnetic scattering rate. Of 
course, large enough values of l/r s ~ h,T c destroy the 
peak and we found the characteristic value of 1/t s , un- 
der which the peak completely washes out, to be ap- 
proximately equal to 0.45T C . The magnetic scattering is 
typically arises due to the magnetic inhomogeneity, re- 
lated above all to A^i-rich clusters, and can be of order of 
the average exchange field in the ferromagnetic alloys like 
Cu\- x Ni x (see I22I and references therein). However, one 
can believe that the magnetic inhomogeneity could be 
weaker in the more diluted alloys with very low exchange 
fields of order of A resulting in small values of magnetic 
scattering parameter l/r s , even compared with T c . 

The appearance of this additional maximum in the in- 
termediate regime etk ~ A is caused by the proximity 
effect in the ferromagnetic region and can be qualitatively 
understood as follows. It is well known that for a diffusive 
SNS junction the zero-bias density of states (DOS) has 
a minigap in the normal region. For low-transparency 
junctions it is A g oc Getk- If the normal metal is re- 
placed by a ferromagnet, the exchange field h shifts the 
densities of states for the two spin subbands in the op- 
posite directions, therefore the minigap in the spectrum 
closes at h ~ A^, but the sharp onsets in the DOS for 
the both subbands survive until they merge the edge of 
the continuous spectrum. So for small enough values of 
A g < h and eV ~ A + h + A g an electron belonging 
to the one spin subband travels through the interlayer 
from one sharp onset of the DOS at —A to another one 
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at h ± A g , while an electron from the other spin sub- 
band can move from —hi A g to the edge of continuous 
spectrum at +A. According to this qualitative consider- 
ation the maximum should be split and the splitting is 
roughly proportional to the value of the equilibrium mini- 
gap A g ~ Getk- We do not observe such splitting at our 
curves. In addition, due to the multiple Andreev reflec- 
tion processes there could arise the subsequent harmonics 
at eV ~ (A + h ± A g )/n. Of course, the proposed quali- 
tative explanation is very crude because it deals with the 
equilibrium minigap instead of real position-dependent 
pseudogap at finite bias and does not take into account 
the behavior of the distribution function in the interlayer, 
which is highly nonequilibrium and obviously influences 
the I-V characteristics. 
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FIG. 2: (a) Possible dependencies of the weak exchange field 
on temperature. The black solid line represents the BCS- 
type behavior of A(T)/T C . (b) The appropriate temperature 
dependencies of the deficit current, normalized to its value 
at h = 0. The inset shows the deficit current eIR/T c as a 
function of the exchange field at T = 0, G — 0.1, eth = 2T C . 

Let us now discuss the behavior of the high voltage 
eV ^> A deficit current Ig e f = I — V/R < 0, where 
R = Rf + 2R g is the full resistance of the junction. The 
inset in Figure [5] represents the typical dependence of the 
deficit current on the exchange field. The deficit current 



is a non-monotonous function of the exchange field and 
reaches a maximum at h ~ A. The non-monotonous be- 
havior of the deficit current leads to very peculiar depen- 
dencies of this current on temperature. The correspond- 
ing curves have the characteristic features determined by 
the temperature dependence of the exchange field. Fig- 
ure dja) shows several examples of possible behavior of 
exchange field on temperature in comparison with the 
temperature dependence of A(T). Figure [2jb) repre- 
sents the corresponding dependencies of the deficit cur- 
rent on temperature. They are normalized to the value 
of the deficit current at zero exchange field. It is seen 
that (i) if h(T) < A(T) in the whole temperature region 

< T < T c , the curve J«fa/(T) = Idef(T)/I def ,h=o( T ) re- 
sembles the appropriate dependence h(T) up to the Curie 
temperature of the ferromagnet Tq u and then follows the 
curve for h — (dotted curve), (ii) in the opposite case of 
h(T) > A(T) for < T < T c I de f(T) monotonously in- 
creases up to T = T c (dashed curve), (iii) if h(0) > A(0) 
and Tcu < T c , then Id e f(T) has a minimum approxi- 
mately at the intersection point of the curves h(T) and 
A(T). For the temperatures higher then Tc u it follows 
the curve for h = (dashed-dotted curve) and (iv) if 
h(0) < A(0) and T Cu > T c , then Idef(T) grows starting 
from the values less then that one for h = 0, intersects 
this line and rises very sharply at T — > T c (gray solid 
curve). 

In summary, we have developed the theory of the I- 
V characteristics of diffusive SFS junctions with highly 
enough resistive SF interfaces. It is found that weak ex- 
change field of the ferromagnet manifests itself only in 
the coherent regime Etk ~ A. For h < A it results in 
the appearance of the additional well-pronounced max- 
imum in the differential conductance at eV ~ (A + h) 
if the magnetic scattering rate in the ferromagnetic re- 
gion is small compared to T c . We have also studied the 
deficit current at high voltages and found that it be- 
haves non-monotonously as a function of the exchange 
field exhibiting a maximum at h ~ A. This gives rise to 
the characteristic features in the temperature behavior 
of the deficit current, determined by the temperature de- 
pendence of the exchange field. We believe, that all the 
discussed features can be measured in the diffusive SFS 
junctions, realized on the basis of weak ferromagnetic al- 
loys. 
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